Photon dressed electronic states in topological insulators: Tunneling and conductance 
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The surface bound electronic states of three-dimensional topological insulators, as well as the 
edge states in two-dimensional topological insulators, are investigated in the presence of a circularly 
polarized light. The strong coupling between electrons and photons is found t o give rise to an 
energy gap as well as a unique energy dispersion of the dressed states, different from both graphene 
and conventional two-dimensional electron gas (2DEG). The effects of electron-photon interaction, 
barrier height and width on the electron tunneling through a p — n junction and on the ballistic 
conductance in topological insulators are demonstrated by numerical calculations. A critical energy 
for an incident electron to tunnel perfectly through a barrier is predicted, where electrons behave as 
either massless Dirac-like or massive Schrodinger-like particles above or below this threshold value. 
Additionally, these effects are compared with those in zigzag graphene nanoribbons and a 2DEG. 
Both the similarities and the differences are demonstrated and explained. 

I. INTRODUCTION 

The unusual energy band structure of topological insulators (TI), as a novel class of quantum spin materials, has 
received a considerable amount of theoretical attention in the last few years—. The energy dispersion is characterized 
by an insulating gap in the three-dimensional (3D) bulk states as well as by topologically protected conducting states 
localized either around the two-dimensional (2D) surface for 3DTIs or around the edge for 2DTIs— . In this paper, we 
adopt the conventional classification (2D/3D) for TIs based on their geometry. 

Quantum spin Hall (QSH) topological states were discovered in HgTe/CdTe quantum wells (QWs). The existence 
of these QSH states is determined by the QW thickness greater than a critical value. For films thicker than 6nm, 
such QWs are exemplary 2DTIs but become conventional insulators otherwise—. Typical examples of 3DTI include 
half-space Bii-^Sbz alloys as well as Bi2Se3, Bi2Te3 and Sb2Te3 binary crystals. The surfaces of these 3DTIs support 
spin-polarized Dirac cones analogous to graphene— 

It has been shown that topological states may acquire an energy gap— i Since it is usually produced as a geometrical 
gap, it requires a finite size along a given direction. The energy gap depends on either the ribbon width for a QW 
in 2DTIs or the separation between two surfaces of 3DTIs. For example, a ribbon width of 200 nm produces a gap 
which is about 0.5 meV. In this paper, we propose another approach for generating an energy gap by coupling the 
topological surface states to circularly polarized light. This type of dynamic gap was predicted in graphene based on 
both a semiclassical approach— and a quantum mechanical formalism— ~—. 

The creation of an energy gap may lead to a metal-insulator phase transition. The conical dispersion of metallic 
graphene has provided unimpeded electron tunneling through a p-n junctions (Klein paradox)— On the other hand, 
Klein-like total reflection has been predicted for bilayer metallic graphene with its massive but still chiral electrons—. 
Chirality is shown to be the key property for total reflection. However, perfect tunneling may still be expected for 
certain values of the longitudinal momentum of electrons in the barrier region because these transmission resonances 
are not affected by the chirality. 

The light-induced energy gap is able to break the chirality— and suppresses the Klein effect in graphene. In this 
paper, we would expect that a similar effect may occur in TIs because the helicity of the topological states is also 
broken by an energy gap. Our numerical results demonstrate a cross-over behavior from Klein-like tunneling in a TI 
to tunneling of conventional a two-dimensional electron gas (2DEG). Here, by Klein- like we mean that the energy 
dispersion of topological states deviates from the Dirac cone. For instance, in a 3DTI, there exists an inherent mass 
term in the effective surface Hamiltonian, which affects the Klein effect in TIs. The interplay between the induced 
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and inherent mass terms, as well as their competing effects on the electron transmission, are the main subjects of our 
investigation. In Ref. [14| , tunneling and transport problems in the presence of tilted uniform magnetic and electric 
fields were studied, and we will briefly discuss the effect of the dynamic gap for these cases. 

The rest of our presentation is organized as follows. In Scc.|TTl we first discuss the dressed topological states and 
obtain their energy dispersion and wave functions for surface states and present an effective surface Hamiltonian as 
well. Similar to graphene, the quantum field formalism predicts a dynamic gap due to electron-photon coupling. The 
distorted valence-band dispersion of the TI is calculated and compared to that for graphene with the main focus on 
broken chiral/helical symmetry. In Sees. IIIII and IIV1 we explore the effect of an induced gap on electron transmission 
through a barrier in a TI for incident energies either less or greater than the barrier height. Specifically, Sec. IIIII is 
devoted to tunneling in 3DTIs along with comparisons to graphene and 2DEG, whereas Sec. lIVI deals with electron 
transmission in 2DTIs to complement the results on electron tunneling in zigzag graphene nanoribbons (ZNRs). 

II. ELECTRON-PHOTON INTERACTION AND DRESSED STATES 

In this section, by including electron-photon coupling, we derive an effective Hamiltonian for surface states of TIs 
based on quantum field theory. Both the single-mode and double-mode optical fields are considered and their energy 
dispersions for dressed electron states are compared. Analogous with graphene-likc masslcss particles, the effect of 
massive particles in TIs on electron states and tunneling are studies. 

Let us now consider electron-photon interaction on the surface of a 3DTI. We obtain the dressed electronic states 
analytically and investigate the tunneling properties of these states. We first assume that the surface of the 3DTI is 
irradiated by circularly polarized light with its quantized vector potential given by 

A = Jo (e + a + e a 1 ") , (1) 

where the left and right circular polarization unit vectors are denoted by e± — (e x ± ie y )/^/2, and e x (e y ) is the 
unit vector in the x (y) direction. The amplitude of the circularly polarized light is related to the photon angular 
frequency ujq by T§ ~ yl/wo- Here, we consider a weak field (energy ~ Tq) compared to the photon energy TtjjJq. 
Additionally, the total number No of photons is fixed for the optical mode represented by Eq. (fTJ), corresponding to 
the case with focused light incident on a portion of an optical lattice modeled by Floquet theory—. 




FIG. 1: (Color online) Transmission amplitude T for Dirac-like surface states for a TI with no quadratic term in the energy 
dispersion, T> = 0. Panels (a) and (b) are density plots of T{4>, e) for a potential V(x) = Vq and width W — 50 nm as well 
as of T(4>, A) for V(x) = X5(x) when e/Vo ~ 1.5, respectively. Panels (c) and (d) show the <j> dependence of T(4>, e) with 
W — 50 nm at e/Vo = 1.5 and e/Vo = 5, respectively, where the red solid curves are for W = 50 nm, the blue dashed curves for 
W — 150 nm, and the green dashed curves for W = 250 nm. 

The non-interacting Hamiltonian of 3DTI was derived in Ref. (l5l ] and we write it as 

n^ = vk% 2] + A^^={fl f k2 ) , ( 2) 
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where I[ 2 ] is a 2 x 2 unit matrix, a stands for the usual Pauli matrices, k = (k x , k y ) is the in-plane surface wave 
vector with respect to the T-point and k± = k x ±ik y . For the 3DTI considered here, the group velocity is of the 
same order of magnitude as graphene, i.e., A ~ Hvf ~ 10 _29 J-m. It was shown that the leading quadratic term 
in Eq. ([2]) is necessary although the higher order terms with respect to 0(k 2 ) may be neglected. The massless form 
of the Hamiltonian in Eq. (j2l with no quadratic term, T> = 0, formally coincides with graphene Dirac cones and 
retains all graphene electronic properties. Specifically, for electron tunneling, these properties include the absence of 
back-scattering for head-on collisions (Klein paradox) as well as distinct tunneling resonances in the electron energy 
distribution. 

The energy dispersion relation associated with Eq. ^ is ef£ l = T>k 2 + (3 A \ k\, where /3 = ±1 is analogous with 
pseudo-spin in graphene. Both A and T> are independent of wave vector k. This dispersion relation shows that the 
particle-hole symmetry is broken by virtue of the massive 2?-term. For completeness, the transmission amplitude 
of the massless topological states (with V = 0) is presented in Fig. [I] Comparing Figs. Ufa) and (b), we clearly see 
a significant difference although the thickness of a potential barrier is only 50 nm. The effect of coupled dressed 
states on the tunneling is much stronger for a (^-function barrier. At the same time, the tunneling resonant peaks are 
broadened significantly compared with graphene. From Figs.[ljc) and (d), we also find that the angular distribution of 
transmission side-peaks at larger angles displays a non-monotonic dependence on the barrier width W for the higher 
scaled electron energy b/Vq = 5. In addition, the broadening of resonant peaks at small angles is also significant in 
comparison with graphene. 

The interaction with the optical mode may be introduced into the Hamiltonian in Eq. ^ via a standard transfor- 
mation of k — > k + eA/h. In Appendix A, we have shown that this transformation leads to the following effective 
Hamiltonian, after the field correction has been neglected, 



H = huj a) a + Vk 2 I [2] + 2QV (k + a + fc_a f ) I [2] + A a ■ k + 2(A (a + a + o _a)) , (3) 

where we have introduced a small parameter £ = eT^j (y2K) to describe the light-matter interaction. We assume 
that the optical mode accommodates a large number No of photons with Wo > 1. Consequently, all the terms ofder 
C 2 ~ O(l/N ) may be neglected. Under these conditions, the energy dispersion associated with Eq. © becomes 

e^' f (A) = N Q hw a + Vk 2 + p^A 2 + (Akf (4) 
with j3 = ±1 and the induced energy gap defined by 



A = ^/W 2 + M 2 - tiwo ~ two ( yj , (5) 

where a = Wo/ (fiuJo) and Wo is the electron-photon interaction energy. For the upper subband with ft = 1 in Eq. (j4|. 
the energy gap is related to the effective mass around k = through 2m* A = ti 2 / [_4 2 /(2A) + T>\ , where the photon 
dressing decreases the effective mass. This is in contrast with single-layer graphene, where electron-photon interaction 
leads to an an effective mass. A similar phenomenon on the effective mass reduction is also found in bilayer graphene 
under the influence of circularly polarized light. The dressed state wave function corresponding to Eq. (U) is given by 

where 7(/3) = Ak/[A + /3y A 2 + (^4fc) 2 ] and <j> = tan~ 1 (k y /k x ). The energy dispersions associated with the Hamilto- 
nian in Eq. ([3]) [also see Eq. (|A17[) ] for two-mode light interaction with electrons are given by 



(A, A) = (n + I) hw + Vk 2 ± yjd(k, A) ± ^£ 2 {K A) , (7) 
where a doubled state space is used for spanning the Hamiltonian, 



Ci(fc, A) = (hujQ/2) 2 + CvAVk 2 + [C, 2 V 2 + A 2 {1 + h/2v 2 )]k 2 , 
C 2 (fc, A) = CM 2 2? 2 fc 4 (I + v 2 ) + Av{C,vV + A)A 3 k 4 
-Zv{2C,V + isA)Ahu Ak 2 + 4(hujo) 2 [A 2 + (Ak) 2 ] , 



(8) 
(9) 
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and v ~ a/2. 




FIG. 2: (Color online) Energy dispersion (red solid curves) for the effective surface model of 3DTI: (a) without light-electron 
interaction [a = 0, T>kp/(hujo) = AkF / (FiuJo) = 0.2]; (b) for single-mode dressed states [a = 0.05 (inner), 0.07 (middle), 0.1 
(outer)] in Eq. Q; and (c) for two-mode dressed states in Eq. Q. Here, the black dashed lines in (a) represent the Dirac cones 
and a parabola, while the black dashed curves in (b) and (c) indicate the asymptotical behaviors for a = 0. 

As QT> — > and uk —> 0, the two-mode dressed states become decoupled and are simply given by 

HN f,Noi}(k, A) = N hLu + Vk 2 ± VA 2 + (Akf , (10) 

e {No+n , No+n} {k, A) = (N + l)hw +Vk 2 ± sj A 2 + (Akf . (11) 

On the other hand, for single mode dressed states, the effect due to the electron-photon interaction is quite similar to 
graphene, except that the energy gap varies as a 2 . However, this dependence becomes negligible under low -intensity 
light illumination. The energy dispersion relations for single and double-mode dressed states of 3DTIs are presented 
in Fig.[2j Comparing Figs.[2ja) and (b), we find that an energy gap is opened at k = due to photon dressing, and 
the Dirac cone is well maintained except for large k values. In contrast, for double-mode dressed states in Fig.jSfc), 
additional mini-gaps appear at the Fermi edge and new saddle points are formed at k = due to strong coupling 
between dressed states with different pseudo-spins. These new mini-gaps and the saddle points prove to have a 
significant effect on electron tunneling. 

The full expression for the 3DTI Hamiltonian which includes the z dependence (perpendicular to the surfaces) may 
also be related to dressing and can be expressed as— (see Eq. (jBjn the Appendix) 

H 3D (k ± , z) = U<i ] D (z) + #$(kj_) . (12) 

where 

*£> W . , C - «!) ?„ + ( <* + ^ f - ^ - {M + M) I + iAA „ ) , (i,, 

uW(k ± ) = -'D ± k% i] -B ± k 2 V 3 ® I[2]+(^ fc ° + o A±k -°^ . (14) 

In this notation, k^ = [k x , k y ), C, Ad, A z , B z , T> z are parameters in the Kane k ■ p model for bulk states, and 
A± , B±, D± are the parameters for surface states. The corresponding energy dispersion relations are linear and 
gaplcss for semi- infinite samples (half-space). However, for a TIs of finite- width, an energy gap is opened due to the 
finite-size effect. We further find that, similar to graphene, electron-photon interaction may modify the energy gap 
due to additional contributions from dressing. 

III. TUNNELING AND BALLISTIC TRANSPORT IN 3DTI 

In this section, we compare surface states of a 3DTI with bilayer graphene in order to find out how w the quadratic 
term in their energy dispersions influences tunneling. The electron tunneling behaves as massless Dirac fermions 
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for high incident energies above a critical value, which has been compared with electrons in a graphene layer. On 
the other hand, the electron tunneling behaves as Schrodinger massive particles at lower incident energies below this 
critical value, which is similar to the 2DEG. Additionally, the modification from a dressed state energy gap to electron 
tunneling is investigated for both low and high incident energies. 

Let us first turn to the tunneling problem associated with the Hamiltonian in Eq. ([2]) in the absence of electron- 
photon coupling. The main input of the p-n junction tunneling is the electron wave function 12 i 16 



* k (x, V) = ^( pt*4> ) e <fc -* +tt »» , (15) 

where f3 — ±1, as before, and <j> = tan -1 (k v / 'k x ) . This wave function is simple and chiral, and an eigenstate of the 
projection of the electron momentum operator along the pseudo-spin direction (chirality/helicity operator) given by 
h = <t • p/(2p). For the tunneling process along the x direction, the transverse momentum hk y is conserved, whereas 
the longitudinal wave vector component k x ^. in the barrier region with V(x) = Vq is determined by 



^=( - ft4+ ^; +4P — )'-"'- « 

Equation (|16[) yields not only propagating but also evanescent modes in the barrier region. Moreover, because of the 
quadratic energy dispersion, both the wave functions and their derivatives need to be matched at boundaries, similar 
to the tunneling problem with respect to bilayer graphene—. 

Electron tunneling and Andreev reflection for the full 3DTI Hamiltonian was studied in Ref.[T3- In the present 
case, we have *4/2? k. Consequently, we may neglect the evanescent contributions. By neglecting the higher-order 
terms, Eq. (fTSf is simplified to 



= / if — Vo ) 2 _ kl . (17) 

' V 22?(e - Vq) + A 2 y V ' 

The result for the special case of masslcss Dirac fcrmions of graphene may be directly obtained from the above equation 
after setting T> 0. 

Clearly, from Eq. (fTTj) . the electron transmission varies with the incoming particle energy as well as the angle of 
incidence. There exists a critical energy e cr <C Vq above which the transmission behaves like Dirac electrons in 
graphene. However, particles with incoming energies below this critical value are transmitted like normal Schrodinger 
electrons in 2DEG. For head-on collisions with k y = 0, the critical energy is calculated to be 



A 2 

Vo ~ ^ • (18) 
22? 



The dimcnsionless two-terminal tunneling conductance, 17(e), may be calculated from 
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ir/2 

s( £ ) = fP = 5 / T(0, £ )cos(4>)#. (19) 

-7T-/2 

The physical meaning of Eq. (|19p is the decrease of the electron conductance in the presence of the barrier due to 
T{4>, e) < 1, where T{(j>, e) is the transmission probability. According to Ref. Go m the case of 2? = may be 

estimated using 



ZlL 

2tt 



dkii 



cosh (k y L x /2) 



(20) 
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where kp is the Fermi wave number and L x and L y are the normalization length and width of the sample. When 
kpL x <C 1, we have tdu\i\v(kpL x /2) ss kpL x /2 and obtain in a straightforward way that Q = (2e 2 /h) (L v /2tt) (ep/A) 
for the Dirac cone with ep — Akp. In the presence of a small energy gap A <C Akp, on the other hand, Qq is modified 
to 



Go= \— J 2^4 ' (21) 



where £f is the Fermi energy at k = kp. 




FIG. 3: (Color online) Two-terminal tunneling conductance for the effective surface model of 3DTI with different barrier widths, 
where X>fcp / A = 0.4. In panels (a) and (b), we present numerical results for tunneling conductance as functions of the incident 
particle energy for the case of Dirac cone and 3DTI. 



The calculated results based on Eq. (Q~9|) are presented in Fig,[3] Comparing Figs.|3][a) and (b), we find that the 
existence of the critical energy in (b) plays a crucial role in the electron tunneling. A series of resonant peaks occur 
above the critical energy e cr in (b) and above the barrier height Vq in (a) although these two values are quite different. 
The former is unique to surface states of a massive particle in 3DTI while the latter is related to the transition from 
Klein tunneling to a regular one for massless fermions in graphene. It is more interesting to notice that the electron 
tunneling below Vq in (a) and below the critical energy in (b) is also qualitatively different, where electrons behave like 
Dirac fermions for Kelin tunneling in (a) and like Schrodinger particle for evanescent wave tunneling in (b). Transition 
when e = e cr is extremely sharp in (b). 

It is very helpful to compare the results obtained in this section with bilayer graphene having quadratic dispersion. 
For bilayer graphene, its lowest energy states are described by the Hamiltonian— 

7i Mg = + 4 , (22) 
2mf, 

where mj, is the effective mass of electrons in the barrier region. In this case, the longitudinal wave vector component 
in the barrier region is given by 



k b 

K x,2 



f3' ^2m b [3(e - V ) - , (23) 



where /3', /3 = ±1. An evanescent wave having a decay rate Kb can coexist with a propagating wave having wave 
vector k x ,2 such that ky + k x 2 = k^ — n\ = 2mb(3(s — Vo). This implies that the evanescent modes should be taken 
into account simultaneously. For bilayer graphene, both wave functions and thei derivatives must be continuous at 
the interfaces. Consequently, the Klein paradox persists in bilayer graphene for chiral but massive particles. However, 
one finds complete reflection, instead of complete transmission, in this case. This effect has direct links with the 
specific electron-hole conjugation, i.e. k x 2 iKj,, in the barrier region—. 

To understand the physics for electron tunneling in bilayer graphene, we first present analytical results for single 
layer graphene, i.e. T> = 0. By taking into account all four modes from Eq. (|16[) . for a (5-potential barrier, we obtain 
the transmission probability as— 

Tg ' 9 = l + sin 2 (A) tan 2 (0) ' (24) 
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reproduces the Klein paradox for the head-on collision corresponding to <j) — 0. However, the periodic dependence of 
the transmission on the scaled barrier strength A in Eq. (|24p is non-trivial. Such perfect tunneling with Ts, g = 1 also 
exists for a set of A values when sin(A) = is satisfied. Furthermore, this prediction is consistent with finite-width 
barrier tunneling [see Fig.Q](c)]. In fact, the expression in Eq. may be derived directly from the general result for 
electron transmission through a very high potential barrier with Vq 3> e, that is, 



COS 2 ((/>) 



1 - cos 2 (k x , 2 W) sin^O) 1 + shr(fc Xj2 W0 tan 2 (0) 



(25) 



where k x ,2 ~ —Vq/A is used for Vq 3> e. It is clear that A in Eq. (j2"4")l plays the role of k x ^W for a finite barrier width 
W. 

For a conventional 2DEG with a (^-function potential barrier—, its trans mission a mplitude is given by 7a = 
2H 2 e/(2h 2 e + m*X 2 W 2 ). Here, we consider a head-on collision with k x = ^2m*e/h 2 and m* is the electron ef- 
fective mass. One may easily see that for A — > (or a very-thin barrier layer) complete transmission (7a —> 1) can be 
obtained. When A — > oo, on the other hand, one gets complete reflection. 
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FIG. 4: (Color online) Density plots of the transmission amplitude T(<j>, e) for bare electron states with energy gap A = in 
the effective surface model for 3DTI with (D/A) k = 0.4. Panels (a), (b), (c) and (d) correspond to fci, ma xVF = 2.5, 1.1, 0.95 
and 0.75. 



Figure [4] presents our numerical results for the transmission amplitude T(4>, e) for several values of barrier width 
in the absence of photon dressing or when the coefficient A = in the energy dispersion relation. From this figure, 
we find that the coupling between dressed states with different pseudo-spins is very strong for a double-mode optical 
field, similar in nature to the result in Fig.[2jc). Here, Dirac-like tunneling above a critical energy may be seen for 
a thick barrier layer, as described by Eq. (|2"5]l . Schrodingcr-like tunneling below the critical energy, on the other 
hand, may only be observed for a relatively thin barrier layer under the normal-incidence condition, as discussed for 
2DEG with A — > 0. With decreased W, electron tunneling below the critical energy is gradually enhanced as — > 0. 
From direct comparison between Fig. [4] and Fig.QJa) we know that the major effect of the massive 2?-tcrm on the 
electron tunneling is associated with the occurrence of a critical energy below which no significant electron tunneling 
is expected for a thick barrier layer. In addition, the critical energy shifts up with decreasing W and the contribution 
from the evanescent mode to the transmission is found finite as long as the barrier width W meets the condition 
W|Im(Ax,a)| < 1. 

We now turn to the model for the 3DTI surface with A / 0. For this, the longitudinal momentum of electron 
dressed states may be approximated by 



Numerical results for the transmission amplitudes based on Eq. (f2l))) are presented in Fig. [SI where the same parameters 
were chosen as those in Fig.|U Comparing Fig.[5ja) with Fig. [Ha) for T> = 0, we find that the effect of a photon- 
induced energy gap ia to produce additional side peaks in the angle distribution of 7~(0, s). After a massive 2?-term 
is introduced, as presented in Figs.[5jb)-(d), the side peaks are significantly suppressed and a critical energy appears. 
At the same time, the Schrodingcr-like electron tunneling below this critical energy is also partially suppressed for a 
thin barrier layer even under the normal- incidence condition by comparing Fig.[5Jd) with Fig.@|d). 
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FIG. 5: (Color online) Density plots of the transmission amplitude T((f>, e) for the electron dressed states in the effective surface 
model of 3DTI with A/Vb = 0.2, T> = (Dirac cone) for panel (a) and (T>/ A) k — 0.15 for all other panels. Panels (a), (b), (c) 
and (d) correspond respectively to fci imax W = 1.5, 1.3, 1.1 and 0.85. 



IV. ELECTRON TUNNELING IN 2DTI AND ZNR 



In this section, we compare edge states of a 2DTI with a zigzag graphene nanoribbon in order to find any similarities 
resulting from an induced energy gap in their energy dispersions. At the same time, the dissimilarities between the 
2DTI and an armchair graphene nanoribbon is also discussed. Additionally, the effect of the decay of electron wave 
functions away from the edges on the electron tunneling are investigated both inside and outside the barrier region. 

According to the Kane model for a HgTe/CdTe semiconductor quantum well— i, the 2DTI system may be effectively 
described by the following 4x4 matrix Hamiltonian, i.e., 



where 



C + M-<p + B)(f4-di) A(k x -id y ) 

A(k x + id y ) C-M-{V-B){k 2 x ~d 2 y ) I ' 



where k x is assumed small. Additionally, we assumed translational symmetry along the x-axis so that k x — > —id x . 
The coefficients A, B, C, T> and M. are expansion parameters within the Kane k ■ p model. Their role in the energy 
dispersion is such that M. is usually referred to as Dirac mass and B is called the Newtonian mass. We note that M. 
changes sign at the critical thinness d cl ~ 6.1 nm. For a layer with d > d c , M. is negative which is referred to as an 
inverted-type TI. 



We note that the Hamiltonian in Eq. (|27[) is block-diagonal, for which we expect two independent wave functions, 
^ + and corresponding to different directions of the electron spin projection. It is straightforward to show that 
\&+ and 'F- are related to each other by a time-reversal operator—, i.e. >F_ = — ia y ^f\. This implies that any 
transmission associated with >F + must accompany another transmission of <F_ in the opposite direction along the 
same edge of the 2DTI. For that reason, the eigenstates <F± are usually referred to as helical edge states—. For all the 
calculations which follow, we made use of the fact that C only appears along the diagonal of the matrix in Eq. (|2"5| . By 
making the replacement e — > e — C, one may show in a straightforward way that the C-term is irrelevant to electron 
transmission. 

Before solving the eigenvalue problem, we estimate— all quantities and coefficients in the Hamiltonian in Eqs. (j2"7| 
and (|28p at d = 5nm (< d CT ). The range of considered wave vectors is chosen as k = 10 _2 A _1 <C kp. In that range, 
the parameters appearing in the Hamiltonian are:— C ~ 10~ 2 eV, Ak ~ Bk 2 ~ T>k 2 ~ \{)~ 2 eV . Consequently, each 
term in the Hamiltonian in Eq. (|2"8"f has the same order of magnitude and none of them should be neglected. Since 
B ~ V, we get B+ = B + V ~ 2B and 6_ = B - V < B. 



Apart from those considerations concerning the quantities and coefficients in the Hamiltonian in Eq. (|27[) , the matrix 
may be separated into two parts. One part depends only on the y— coordinate whereas the other one is fc^-dependent. 
The eigenvalues are determined from the following secular equation 



9 



(M B + {kl - e) -s)(M- B.(k 2 x £ 2 ) - e ) = -A 2 (k 2 x £ 2 ) , (29) 

where £ characterizes the decay rate of edge state electron wave functions with a finite width in the y direction. After 
neglecting the i3_-term, we obtain from Eq. (|29|) 



* * A 2 -B + {e-M) ■ {M} 

The electron wave function, tunneling and transport properties for a system with finite width differ significantly from 
those of a semi-infinite model because of a gap in the energy dispersion. 

In the limiting semi-infinite geometry for a 2DTI, we obtain the following exact ID effective edge model along the 
y-direction,— 



Uf^ = Ak y V z , (31) 

whose corresponding dispersion relations arc e^™ 1 = ±Ak y with corresponding wave functions are the eigenfunctions 
of <r z , i.e., — {0, 1} and = {1, 0}, which yield a transmission amplitude T ee m = 1 through a barrier of any 
height analogous to the Klein paradox for head-on collisions in graphene. In contrast to the 2D model in Eq. (f2"T)l . we 
find that the electron transmission for edge states with a finite width is substantially suppressed. 

Taking the limit L y — > oo, we can also consider bulk states which are located far away from either edge. In this 
case, we make the substitution k y = —id y . By retaining terms up to order C(fc 3 ), calculation leads to a new energy 
dispersion for A 3> Bk 



£■ 



bulk 
2D 



= -Vk 2 ± y/M 2 + (A 2 - 2MB)k 2 . (32) 



The Hamiltonian in Eq. (|2"5|) assumes the the simple form ±A4 tr z at k = 0, where the gap parameter M. strongly 
depends on the thickness d of the quantum well and can be arbitrarily small or even set equal to zero. The energy 
dispersion relations in Eq. (|32p formally reduce to those of gapped graphene or graphene irradiated with circularly 
polarized light. The tunneling problems for this case were addressed in Refs. (l~8l | and [HI]. The most significant 
effect of radiation on electron tunneling is the breaking of chiral symmetry on the order of 0(A4 2 ). Consequently, 
significantly different behavior in the electron transmission at small incident angles (\k y \ -C \k x \) is expected compared 
to infinite graphene in the absence of light illumination—. 

As derived from Refpi, the decay rates £12 introduced in Eq. ([2^]) are given by 



Ay/1 - (V/B) 2 y/A 2 (B-V)-4BM(B + V) 

2{B + V) 2BVB + V ' 1 ' 

The above equation, in conjunction with Eq. (|29[) . leads to the following energy dispersion relations 



MV . 



^{B 2 -V 2 )k x 



o{kl) , 



(34) 



where A z is defined by 7 



AAM(B 2 -V 2 ) 

(35) 



z B 3 [A 2 B-4M(B 2 -V 2 )) ' 
Additionally, the transposed wave function associated with £^2 under the limit of 3> 1 is 

= ^{/c - M, ( g + p )fe + 6) (/o + p 8hU)} (36) 
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with the notation j3 = ±1, Sk = ±1, Af is the normalization factor, and we express the wave function components in 
terms of odd (/ Q with — ) and even (/ c with +) functions, defined by— 

e fiy -|- e -£iv e &v ^ q-&v 

fe(o)(V) = e ^ lLy /2 ±e -^L y /2 ~ e bL y /2± e -&L y /2 ' ( 37 ) 



We now turn to a comparison of the electron tunneling in 2DTI with graphene zigzag nanoribbon (ZNR) quasi-lD 
edge states. In Ref. (22|, the authors calculated the ZNR wave functions as 



Be 



S mh[X v>n (L y /2 + y)} 
sinh[A y , n (L y /2 - y)\ 



(38) 



where X y n is related to the edge distribution of electron wave functions. The boundary conditions are such that each 



of the wave function components vanishes at one of the two ribbon edges (zigzag configuration) . 
following relation 



This leads to the 



k x Ay 



t/.n 



k x ~\~ A* 



(39) 



Equation (j3"9")l demonstrates that the real solutions for X yn exist only if both 
are satisfied. We will assume the ribbon is sufficiently wide so that k x L y 
real \ y , n , the edge states decay exponentially with decay length 1/X y<n . 
2DTI discussed above in this section. However, the traverse wave function 
a graphene armchair nanoribbon (ANR) which has a plane-wave type wave 
from Eq. (f39|) that for large k x one gets X y . n 0, i.e., fast moving electrons 
negligible decay rate. Therefore, the wave function in this case extends far 
ribbon width, there exists a maximum value for the decay rate. 



the conditions k x > X yn and L y X y ^ n > 1 
> 1 for k x in the range of interest. For 
This situation is similar to the case of 
distribution is drastically different from 
function ^ANR(y) ~ e lkyV . We also note 
with large longitudinal momenta have a 
away from the ribbon edges. For chosen 




T=l 
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2 
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FIG. 6: (Color online) Zigzag nanoribbon wave function and tunneling conductance. Panel (a) demonstrates transverse decay 
rates A Sj n depending on the longitudinal momentum k x for various nanoribbon widths. Panel (6) shows the y-dependence of 
the edge state wave functions inside the barrier region (A;,) and outside the barrier region (Ao), as well as contour plots of wave 
functions. Panel (c) presents density plot for the dependence of two-terminal conductance g(e) on Ao and A;,. 



The energy dispersion for ZNR corresponding to Eq. eqrcfznrl is given by e? NI 



\2 



We confine our 



attention to low potential barriers, such that \ Vt n is real within a barrier region. Unimpeded tunneling in ZNR was 
investigated in Ref. [23| ■ Additionally, we calculated the two-terminal tunneling conductance by making use of 



9(e) 



Q{e) 



1 

Li, 



Ly/2 



T(e, y) dy 



(40) 



-Ly/2 



showing how the ribbon conductivity is modified by a barrier region. 
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From Fig.|6ja), wc find that the wave-function transverse decay length A~^ decreases with increasing width of a 
zigzag nanoribbon. Figure [Bfb) show us that the minima of edge-state wave functions in the transverse direction are 
at the center (y = 0) for both inside and outside barrier regions. Furthermore, the wave functions are symmetric with 
respect to y = 0. Our numerical results based on Eq. (|4T))) are presented in Fig. [51c). The diagonal line in Fig.[SJc) 
reflects the fact that T(e, y) = 1 [or g(e) = 1] for the same decay rates A = A;,. The periodic resonant peaks with 
respect to A& can be clearly seen in Fig.[6tc) similar to the prediction by Eq. (f24l . 

V. CONCLUDING REMARKS 

In summary, we have analytically obtained the energy dispersion relations as well as the wave functions of elec- 
tron dressed states in TIs irradiated by circularly polarized light. A number of helical systems, such as graphcnc, 
nanoribbons and topological insulators, are compared . Similar to graphene, the electron-photon coupling in TI leads 
to a energy gap in the electron energy dispersion relations and eigenstates with the broken chirality symmetry. The 
tunneling over a square potential barrier is modified significantly if there exists a gap in the energy dispersion of 
electrons. The combination of Schrodinger-like (massive) and Dirac-like (masslcss) electrons tunneling below and 
above a critical energy gives rise to very novel properties. We have further found that some lower energy subbands 
become nearly dispersiveless with increasing light intensity, which results in unusual electronic properties non-e in 
graphen 

As shown in Ref. Q for graphene, laser power 10 2 W may produce an energy gap A ~ 10 — 100 mcV required for 
making the effect noticeable for THz light frequencies at room temperature. Since the linear-term coefficient (group 
velocities) of the 3DTI Hamiltonian has the same order of magnitude as graphene, we expect that an experimental 
verification to be possible. 

Although the linear term in the energy dispersion of 3DTIs gives spin-polarized Dirac cones with chirality of the 
corresponding eigenstates, the appearance of a quadratic term introduces additional effects on electron tunneling. 
Our calculations for electron tunneling through a square potential barrier indicate that electrons may be transmitted 
either like chiral particles as they do in graphene or like conventional 2DEG electrons, depending on the incident 
particle energy above or below a critical energy. 

We have also investigated tunneling properties of edge states for massive particles in 2DTIs, such as HgTe/CdTe 
quantum wells, where an energy gap is introduced by the finite width of samples—. From the analysis of dressed 
edge-state wave function in the presence of circularly polarized light, we have found that particles can always freely 
propagate along the ribbon edges, which is similar to zigzag graphene nanoribbons but different from semi-infinite 
2DTIs where either perfect transmission or complete reflection is obtained for electron-electron and electron-hole 
transitions, respectively Although the wave functions in both 2DTI and zigzag ribbons are localized around the 
edges, the decay rate of edge state wave functions of 2DTI does not depend on the longitudinal wave number, in 
contrast to the decay rate in a zigzag ribbon. 
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Appendix A: Model of 3DTI Surface Irradiated with Circularly Polarized Light 

Tthe Hamiltonian describing the surface states (at z = 0) of a 3DTI to order of 0(k 2 ) is given by 

where k = (fc^, k v ) is the un-plane surface wave vector and k± = k x ± ik y . The energy dispersion associated with this 
Hamiltonian is given by = T>k 2 + f3Ak with f3 = ±1. 

We now turn to the case when the surface of the 3DTI is irradiated by circularly polarized light with vector potential 
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A = Jo (e+a + e_a t ) , (A2) 

where e± = (e x ± ie y )/\/2, e x and e y are unit vectors in the x and y direction, respectively. Consequently, the 
in-plane components of the vector potential may be expressed as 

A x = ^(a + tf), A v = i^(a-tf). (A3) 
In order to include electron-photon coupling, we make the following substitutions for electron wave vector 

k x — > k x + —— = k x + —7=-{a + a 1 ) , 



ky > ky + ^ = ky + ^^^fl ~ ^) 



V2eJ"o „ t %/2eJ" „ 
fe+ — > fc+ H — a 1 , fc_ — > k- H — a , 

n, ft 

k 2 = k + k_ ^k 2 + ^^{k + a + k^)+ [^P^j ^a. (A4) 

In our investigation, we consider high intensity light with Nq = (aJa) ^> 1, and then, dd^ ~ d^d due to dd''' = d^d + 1 

for bosonic operators. We adopt this simplification only for the second-order terms ~ (v^eJ-o/d) but not for the 
principal ones containing Hluq. With the aid of these substitutions, the Dirac-like contribution to the Hamiltonian in 
Eq. (|A1[) becomes 

%Dirac = ' k = -A 0?-k+ + CT + fc-) 

^A(V-k + + V+k-) + ^?^A(V-tf + V + a) , (A5) 

n 

where <?± = (ct^ ± ia y )/2. To describe a full electron-photon coupled system, we also need to take into account the 
photon energy term huio d^d. This yields 

H = (huj a + 4£>C 2 ) a f d + Vk 2 I [2] + 2(1? (fc+d + k-a)) l [2] 

+A(a+k- + a-k+) +2£A(a + a+a^) , (A6) 

where £ = eT^ji^piK). We may also rewrite the Hamiltonian in Eq. (|A6[) in matrix form as 

H = (fiujQ + 4£>C 2 ) ^ a + 1 + 2 + 3 

M (4 -^ +4+ V/ + ^)) +2 K°>o). < a7 > 

where i = 'H^ i denotes the initial surface Hamiltonian with no electron-photon interaction, 3 gives the principal 
effect due to light coupled to electrons (the only non-zero term at k = 0) and 2 is the leading term demonstrating 
the difference between dressed states in graphene and 3DTI. 

We know from Eq. (|A7[) that the Hamiltonian at k = reduces to the exactly solvable Jayness-Cummings model, 
after we neglect the field correction on the order of C(C 2 )- We obtain 

"H k =o = dw d t d + 2C-4(a + d + v-a)) . (A8) 
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Following the method adopted in Refs. [9J and [13J, we expand the eigenfunctions of Eq. (|A8|) as 

l*?,JV„> = fH,No I t, N ) + Vf,N \l,N + l) , 

\Kn ) = Vl,No I I, ^o) - »i,n 1 t, N - 1) . (A9) 

By using the properties 

at|tl,iVo) = v^VoTT|t4,iVo + i) , 

a|n,^Vo) = v / ^oin,^o-l> , 
?± I It^o) = (1 - «t,+)(l - <^,-) I U,JVo> , (A10) 

we obtain the energy eigenvalues 



4 i 

1 /l 1 





" 2 /\r 


1 










1) 




-2± 





iVo± 2 T (i + i a j = NoT ~ ' (A11) 

where a 2 = 2(AN /(haj ) with AT S> 1. The energy gap at k = has been calculated as A = e°_ — e+ « (a 2 /2) Sw - 
We note that there is no difference between graphene and the surface states of 3DTI at k = 0, and therefore, the 
result in Refs. [j| and [24| are relevant to each other. 



The expansion coefficients in Eq. (|A9j) arc calculated as 



Mn,w - J ov n,No , v n,x - \ ov n ' N ° ) ( A12 ) 





i 








4) 


" 2 





where 



^ = yi-^-^ + - ± -j. (A13) 

In all further calculations in this Appendix, we assume a <C 1 and iVo 3> 1, corresponding to a larger number of lase 
photons but weak light coupling to electrons as a perturbation to the electron energy. This leads to the following 
approximate expressions 

a 2 

m,N — fJ-i.No — m,N +i — Vl.No+l — 1 , (A14) 

^t,jv ^ n.^o — ^t.JVo+i — n,N +i — 2 ■ ( A15 ) 

Consequently, we only need to keep one pair of index- free coefficients {/i, v\ such that fi = cos($), v = sin(<£>) with 
$ = tan -1 (a/2). Furthermore, from the first equity in Eq. (| Al 1 j) . it follows that the energy gap still depends on JVo 
in general. However, the difference A% +1 — A^ ~ A /No is so small that we neglect it for Nq 3> 1. 

Generalizing Eq. (|A9[) . we still expand the wave function over the eigenstates of the Hamiltonian in Eq. (|A8[) for 
k ^ 0, i.e. 

JVo+l 
j=N t4 



*?,iVo>+ S Ui*liVo>+ 3 U+ll*U+l)+ S U+li^ ) iVo + l) ■ (A16) 
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Using Eq. (|A16[) we project the full Hamiltonian in Eq. (IA6|) onto the representation {3^ JVo ,S^j Vo ,H^j Vo _ | _ 1 ,H^j Vo+1 }. 
This yields 



/ N hu + Vk 2 
H 2 Ak+ 



A n 2 Ak- 
N huj + Vk 2 



vfiAk- - 2(Vk_ 
V -2ufj,Ak- - 2QVk^ 



\ivAk+ - 2QVk+ 




\ 

-2v(iAk+ - 2QDk + 



{N + l)hu + Vk 2 
[i 2 Ak + 



[i 2 Ak- 
{No + 1) hwQ + Vk 2 + A / 



where A « A is the energy gap at k 7^ and we have employed the relations for N = Nq or Nq + 1 

£ ■ k |*^) = £ ■ k t,N) + v\l,N+l)} = fdt+\ i, N) +vk_\ t,iV+l) , 
V ■ k|*J w ) = 9 ■ k ( M | |,iV> - v\ t, JV - 1)) = nk-\ t, AO - vk+\ i,N-l) , 

= y/N + l\* n , N+1 ) . 



a f \*°n, N ) 



(A17) 



(Ai8) 

(A19) 
(A20) 



After calculating the eigenvalues for the Hamiltonian in Eq. ()A17j) . we obtain closed form analytic expressions with 
energy dispersion 



(k, A) = (n + ^jhoj + Vk 2 ± y/d(k, A)± VC 2 (fc, 

Ci(k, A) = (fi^o/2) 2 + (vAVk 2 + [( 2 V 2 +A 2 (1 + 5/2v 2 )] k 2 , 
C 2 (fc, A) = ( 2 A 2 V 2 k 4 (l + v 2 ) + Av{tvV + A)A 3 k 4 
-iv{2C,V + vA)Atuv Ak 2 + 4?i 2 w 2 (A 2 + A 2 k 2 ) . 

Taking the limits QD —¥ and vk — > 0, we are left with two uncoupled energy subbands 



A) , 



(A21) 
(A22) 

(A23) 



£ {Not,N i}(k, A) = Nohuo + Vk 2 ± ^/A 2 +A 2 k 2 , 



e{*b+it,JVo+U}(*> A) = (JV + 1) tiu + Vk 2 ± ^/ A 2 + A 2 k 2 . 



(A24) 
(A25) 



Therefore, we conclude that the effect of electron-photon interaction is quite similar to graphene as far as one photon 
number Nq is concerned. The main difference being that the energy gap in the 3DTI is of order 0(£ 2 ), which may be 
neglected for low intensity light. Consequently, the energy dispersion relation becomes 



e p (k, Ac) = N hco + Vk 2 + PyJ [A + 0{C 2 )Y + (Ak) 
where (3 = ±1 and A is the photon- induced energy gap as in graphene. 



(A26) 



Appendix B: Dressed 3DTI Electron surface states 

It follows from Appendix [5] that the effect due to ~ A 2 -terms would play a the role only if two optical modes are 
considered. This effect is of order 0(( 2 ) and may be neglected in our calculations. When the surface of a 3DTI is 
irradiated with circularly polarized light, the light will penetrate into the sample and decay exponentially away from 
the surface, similar to the wave function for a surface electronic state. Therefore, by including this decay effect, the 
vector potential in Eq. (|A2[) is generalized as 

A = T (e+a + e_a f ) e"^ , (Bl) 
where l/£ is the decay length. By including the z dependence, the Hamiltonian of the electron system is found to be 



(B2) 
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where 

ti&iz) = (C- V z dl) I [4] +{ {M+ B ^ ^ iAA ^ ° „ ) , (B3) 

H^±) = -V±#lw-B±k>V 3 ®I m +( . ° . (B4) 

\ .Aj_K-|- (Tl U / 

Here, = (fc K , k y ), C, A4, A z , B z , T> z are parameters in the Kane k • p model for bulk states, and A±_ , B±, T>± are 

parameters for surface states. We note that the parameter C may be eliminated by simply shifting the band edges. 

Since the electron-photon interaction occurs in the region close to the surface, we only need to consider (kj_) in 



Eq. (|B2j) . which contains the coupling with the incident circularly polarized light. 

When k = {kj_, k z } = 0, the Hamiltonian in Eq. (|B4|) plus the single-photon energy as well as the light-electron 
coupling together give 

a 

■H k =o = foo a^a + 2(A± I t a I . (B5) 

o at Q 

a< 

By introducing a small dimensionless parameter a' = 2(Aj_/M <C 1, and the two matrices 



1 
1 o 




r+ = o o ^ | , r_ | _ u u | , (B6) 



the Hamiltonian in Eq. (|B5|) may be rewritten compactly as 

W k =o = 7^oa t a + a / (r-a + + r + a) . (B7) 

In analogy with Eq. ((A9|, we construct a basis set containing four states, i.e., | + + ) = |1, 0, 0, 0), | H — ) = |0, 1, 0, 0), 
| — h) = |0, 0,1,0), | ) = |0, 0, 0, 1). For this basis set, it is a simple matter to show the following properties 

a) \N, ±±) = y/N + T \N + 1, ±±), 
a\N,±±) = VN\N- 1,±±), 

r + |iv,+±) = o, r+|Ar,+±) = o, 

f + \N,-±) = \N,+ T ) , T_\N,+ T ) = \N,-±) . (B8) 

We may also expand the dressed electronic states at k = over this basis set leading to 

= W\r,+I + +,N) + p N , + \ + -,N) + u N+1 ^\-+,N + l) + v N+1 -\ --7V + 1) , 
|*+-,jv) = ^,+ | + +,iV) +mat,+ | + -,N) + u N+1 ,-\ - +,N+l) + u N+1 -\ - -,N+1) , 
\^-+,n) = "n,+ \ + +,N) + p n ,+ \ + -,N) + Hn+i,-\ -+,N+1) + i/at+1,-1 - -, N + 1) , 
\y a —, N ) = vn,+\ + +,N)+ v N ,+\ +-,N) + i/jv+i.-I -+,N + 1) + Hn+i,-\ - -,N + 1) . 
For N 1 and a' < 1, calculation shows that 

I \ 2 

a 



^n,+ — fJ-N- — MA r +i,+ — Mw+i - — 1 — ^— ) , (B9) 

a' 

^ fjv,- - VN+1.+ ^ ^w+i.- ^ — • (B10) 
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Introducing the pair of operators 



Tl2 = 



we obtain 



ria|JV,- 

r 34 \n, - 







o 1 
o o 




\N, 



r 34 = 



o o 

1 











r 34 |iv,++) = \n,+-) 

T 12 \N,—) = \N,-+) 



For the Hamiltonian in Eq. (|A7I) within the 2x2 subspacc, as a special case, we introduce 



(Bll) 



(B12) 
(B13) 



c+ = Ti2 , cr_ = r 34 . 

Finally, we can rewrite the full Hamiltonian in Eq. (|B2[) for k ^ 0, using the basis set \Nq, ±±), as 



(B14) 



"H 3 D(k^, z) = N huj - 



C+ - V*,-dl - T>± 


_fc 2 - A 




-iA z 


d z 





A±k- 


— iA z d z 


C- 


- v z 


+ Z ~ 


"D-L. + k 2 - A 


A±k- 





A ± k + 


c+ 




-91- 


T> ± _k 2 + A 


iA z d z 




A±k + 











iA z d z C- 


-V z , + d* + V ±: + k 2 + A 



(B15) 



where the ± signs correspond to opposite pseudo-spins in the basis set. Making use of Eq. (|B15[) . we obtain the secular 
equation for the energy dispersion relations when A = 0, that is, 



s 2 - A\k 2 + NiNi + A 2 z i 2 + £>_Ni£ 2 + V + n 2 £ 2 +V_+ V + N 2 £ 2 + V + V^ A - (Ni + N 2 + 2£>£ 2 ) e = , (B16) 
where we have introduced the notations 



Ni, 2 (fe) =C±M + (V ± TB^)k 2 , 
V± = V z ± B z . 



(B17) 
(B18) 
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